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Abstract 

We start with comparisons of hierarchies in Biology and relate it to Quan- 
tum Field Theories. Thereby we discover many similarities and translate them 
into rich mathematical correspondences. The basic connection goes via scale 
transformations and Tropical geometry. 

One of our core observations is that Cellular Automata can be naturally 
introduced in many physical models and refer to a (generalised) mesoscopic 
scale, as we point it out for the fundamental relation with Gauge Field Theories. 

To illustrate our framework further, we apply it to the Witten-Kontsevich 
model and the Miller-Morita-Mumford classes, which in our picture arise from 
a dynamical system. 
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1 Hierarchies in physics and biology and struc- 
tural similarities 



One of the most intriguing questions is, how can simple constituents generate complex 
structures or patterns and how are they mathematically described. 

In the physical domain a prominent example is given by particles, e.g. atoms, which 
serve the role of the basic building blocks, to yield compound objects, e.g. solids or 
liquids. 

Usually one does not only make the distinction between basic constituents and com- 
plex compound systems, but one also associate characteristic scales (e.g. energy, 
length, time) to the objects themselves where they arise. In the example of atoms, 
one speaks about the micro as compared to the macro world. Guided by experience, 
the mathematics which is capable to describe the objects at the different scales, also 
changed over time. Probably the most prominent example in this direction is the 
transition from classical to quantum mechanics. 

However, one also considers an intermediate range, usually termed mesoscopic, 
which is situated "in-between", both from the metric point of view but also from 
the mathematical description. 

But, as mesoscopic structures build a bridge between microscopic and macroscopic 
objects, mesoscopic systems inherit characteristics of both worlds, e.g. classical and 
quantum signatures. 

Mathematically, in the situation of quantum mechanics and classical mechanics, a 
"semi-classical" description or a "WKB approximation" can be applied. Other pos- 
sible ways to describe this transitional regime are the path integral formalism, with 
a scaling parameter, e.g. h, but also random matrices. 

Molecular biology is considered to be a fundamental theory to analyse many of the 
biological phenomena occurring at, e.g. the cellular level. The method itself is based 
on interaction systems among a set of words, build on an alphabet of four letters, the 
DNA, and the mathematics involved is mainly of an algorithmic and combinatorial 
nature. 

A particular aim is to understand how from such a mechanism macroscopic patterns 
arise, as e.g. the stripes of the skin of fishes. 

In principle, these should be explicable from a physical point of view by molecular 
interaction systems and their electrons, but nevertheless it is mysterious that even 
though the microscopic behaviour appears to be random, after several interaction 
steps at various time scales, macroscopic patterns emerge. Henceforth, we shall refer 
to such systems of dynamical structures as the dynamical topology of scale transfor- 
mations. 

From a dynamical point of view, there are three distinguishable hierarchical steps. 
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micro, mesoscopic and macro. However, compared with the usual physical situation, 
as described at the beginning, there are differences concerning the mathematical 
aspects of the dynamics. Whereas in physics differential equations are used as basic 
mathematical tools, for biology the automaton is fundamental. 

However, what is common to both is that one finds a relative structure, i.e., the 
micro- mesoscopic- macro hierarchies. 

As we have already pointed out, a structural similarity exists between biological 
and physical systems. 

One of the common features of their respective systems is that both consist of hier- 
archies of systems of different sizes. The sizes of these hierarchies of three types in 
biology and physics are mutually very different. In many interesting cases the sizes 
of the biological hierarchies are much bigger than those of the physical one. So the 
intrinsic nature of the individual materials will be different from each other. What is 
similar, is the global structure of these systems, i.e., the mechanisms of creating such 
hierarchies. Such basic hierarchies in both biology and physics, contain many refined 
ones. In order to compare such biological and physical hierarchies systematically, it 
will be natural to categorise these structures, with the similarities interpreted as the 
existence of specific functors between them. So, a particularly important goal for us 
is to construct a category of biological hierarchies, but also the associated functor 
between the physical and biological realm: 

Category of biological hierarchies Category of physical hierarchies. 

One of the conclusions from our claim of the existence of a functor between biological 
and physical hierarchies is that many physical systems should posses reasonable dis- 
crete, or characteristics of a computer (cf. [5J). This is the beginning of our program 
to discretise several physical systems, guided by the comparison with the biological 
side. 

So, in later sections, we shall study some concrete comparisons, based on mathemat- 
ical realisations, coming from both physics and biology. 

2 Mesoscopic models 

We shall describe now two models which are mesoscopic by their nature, however one 
of them is usually not primarily seen from this perspective. But looked at it this way 
one will immediately be rewarded mathematically. 

(A) Gauge Field Theories: A category of physical objects that naturally fits into 
our general framework are Gauge Field Theories and instantons. 

As it is known [191 [IB] , gauge theories can be discretised yielding lattice gauge field 
theories, which, as e.g. at the sub- nuclear level in Mesoscopic QCD {[i3\) posses also 
a mesoscopic regime. 
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Further, the instantons corresponds to field configurations that locally minimise the 
classical action, with the action functional decomposing the path space it into topo- 
logical sectors, corresponding to topological charges. 

Now, in the semi-classical description, the wave function can be expressed by use of the 
Maslov indices of its trajectories under the constrains on the Lagrangian submanifolds 
([!]). It is understood as a Morse theory on the path spaces Q{M,p,q) which 
are infinite dimensional over Riemannian manifolds. By taking the energy functional 
as a Morse function, one can trace a parallel theory to the finite dimensional case. 
The Morse indices tell us the structure of the homotopy types of the path spaces 
Q{M,p,q), and it gives information on the topology of the underlying manifold M 
(see PI). 

Let M be a closed manifold and / : M — )• M be a Morse function on it. To each critical 
point, the degree of the Hessian is assigned, and it is called the Morse index. In fact 
it is a special case of the Maslov index. Now, Witten's complex (C(/)*, d) is a chain 
complex whose chains of degree k are generated by the sets of critical points of Morse 
indices k. The connecting orbits are given by the gradient fiow of /, which yields 
the boundary operators for the complex. Its homology is canonically isomorphic to 
the ordinary homology H^{M : Z). This is a striking point where Morse functions 
connect dynamical structures with topology, passing through the Maslov indices. 

One has a loop space version of this, which is called symplectic Floer homology [1]. 

Let (M, a;) be a symplectic manifold, and choose a compatible almost complex struc- 
ture J on {M,u). On the loop space C°°{S^ , M), one obtains an action functional 
(pH '■ C°°{S^ , M) — )■ M when 7T2{M) = 0, by use of a Hamiltonian function. Now, 
the action functional (pn plays the same role as the Morse function above, and so 
one obtains Floer's chain complex {C{M, H)^,d) on the loop space, which admits 
richer structures than just homology groups. The critical points of the action are 
the periodic orbits, and for each periodic orbit, one can associate the degree by the 
Maslov index. Thus as the case of Morse functions on manifolds, one can construct 
the Witten complex, and obtain the Floer homology HF{M, H, J). When one for- 
gets the extra structure mentioned above and regards its homology HF^:{M, H, J) as 
abelian groups, then it is canonically isomorphic to the singular homology Hi,{M\ Z2) 
of M. Thus one obtains semi-classical and classical objects from the forgetful functor: 
HF,{M,H, J)^H,{M; Z^). 

(B) Iteration dynamics for families: In order to understand the mechanisms of 
biological systems, we propose to express such phenomena as iterations of families of 
maps, and to interpret them as microscopic orbitals. So, let us take two interval maps 
/o, fi : [0, 1] — 7- [0, 1], and denote their iterations as Of-{x) = {fi'i^)}'n=o^ i = 0,1. 
We may regard them as oscillations of a molecule, and therefore we can formulate 
the interaction between Of^ and Oj^ in the following way: let X2 = {{ko, ki, . . .) : 
ki e {0, 1}} be the one-sided full shift with alphabet {0, 1}, and choose any element 
k = {ko,ki,...) G X2. Then, with respect to k E X2, vfe define the interaction 
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between Of^ and Of-^ as a family of interval maps {/i" : [0, 1] — )■ [0, l]}^o where- 
h\x) = fk, o fk^_^ o ■ ■ ■ o fkoix), 1 = 0,1,... . 

Although the family of maps {/i"} happens to behave in a complicated manner, by 
discarding some of the information via a projection map, patterns can be observed 
to emerge. So, let n : [0, — i- {0, 1} be the projection defined as 7r([0, |)) = and 
7r((|, 1]) = 1 and then consider the assignments for a.e. x G [0, 1]: 

k' = 7c{h%x), h\x), ...)) = (tt o h^{x), 7foh\x),...)eX2. 

Thus for each k, one associates another element k' G X2, which we call the interac- 
tion map: 

^x) = ^x,fo,fi):X2^X2. 

Their time scales are much larger compared to the n of the oscillation {h^}, as they 
can be determined only after knowing all the values n o h^[x) up to n = 00. So, 
according to our definition, this is a mesoscopic scale. 

Now, if we write $(x)(A;o, ki, . . .) = (A;q, k[, . . .), for fcj, k^ G {0, 1}, then although k'j is 
in general determined from the data ko, ki, . . . , kj, in some cases it is already given by 
knowing only the data kj-m, kj-m+i, . . . ,kj for some m independent of j. This prop- 
erty turns out to be quite unstable and will be broken even for small perturbations of 
the original maps /o, /i. As such a dynamics corresponds to an automaton, the orig- 
inal interaction can be regarded as a deformation of automata, whose mathematical 
structure is essentially discrete and finite. 

As next we will show the relation of $ with integrable systems. In fact, to a variant of 
the above interaction maps and interval maps, the associated flow {$(x)*(A;)}^o yields 
a solution of the Lotka-Volterra cell automaton, which arises from a new method to 
obtain a scaling limit, namely Tropical geometry. This geometrical theory permits to 
reduce discrete dynamics to cellular automata, passing through specific scaling limits. 

Specifically, Tropical geometry associates rational functions parametrised by t to 
(max, +)-functions, which are piece-wise linear functions on and it transforms 
the dynamics of the original PL function into a complex, parametrised dynamics on 
the affine algebraic varieties Vt C C^. 

As the scaling parameter tends to infinity, the phase spaces of the dynamical systems 
have to be extended to CP^ D C^. At t = 00, the dynamics on CP^"^ = CP^\C^ 
looses the detailed information it had for finite values of t. 

Given the fact that Tropical geometry connects integrable systems, described by 
PDEs, with cellular automata, our perturbation via interactions of maps will allow 
us to extend such a connection to 'neighbourhoods' of the dynamics, including non- 
integrable systems pT] . 

Now, for the hierarchies of scaling limits, i.e.. 

Micro =^ Mesoscopic =^ Macro 

{/o ^ {/r(2/)}m =^ Automaton =^ Solitons , 
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the interaction given by a piecewise linear map on R^, reduces to a cellular automaton 
{A : Vi + max(0, V2 + V3) = V2 + max(0, + V4)}, which can be transformed, again by 
Tropical geometry, into a rational dynamics on the three dimensional affine algebraic 
variety V = {( 21, 22, -2^3, -2^4) : -^2 + 2^1-22-24 = Zi + Z1Z2Z3} C C. Then, as it is shown 
in |10], one can transform the above dynamics into a KdV flow, by applying Hirota's 
classical method |s8J for going from the discrete KdV to the continuous KdV equation. 

Now, let us discuss mesocopic aspects of Tropical geometry in a wider context. 

A particular phenomenon in biology is that among scaling parameters from finite 
values to infinity, many biological systems work 'quite' stably on an intermediate 
range between small values and 00, and even under changes of it. But at t = 00, 
such systems will change very differently. For temperature many biological systems 
will break down at t = 00. For hierarchies of biological systems, scaling limits create 
macro patterns given by very different dynamics from molecular interactions. 

Now, let (/} be a piecewise linear map given by a (max, +)-function. Then Tropical ge- 
ometry accociates with it parametrised rational maps ft, t G (1, 00), which are related 
to (fhy ft = \og^^ oLptoLog^, where the (ft are approximations of for limt^oo = ^■ 
The correspondence between (f and ft is one-to-one for their respective presentations. 

We shall now give the following examples as illustrations: (1) Let us consider the 
(max, +)-function if{x, y) := max(max(0, y) — x, —x) which corresponds to a rational 
function g{z,w) := and which happens to be time-independent. 

The dymanics given by x„ = (f{xn-2, Xn-i) is for any initial value recursive with 
period 5. Thus one might expect that another dymanics, Zn = g{zn-2, Zn-i), will be 
the same, since two dynamics given by g and (ft are mutually conjugate and the limit 
limt^oo V^t = holds. However as a straight-forward calculation shows, this is not the 
case [To] . 

Mathematically, such phenomena occur, because the limit t = 00 does not exist for 
the conjugate maps Logt, and so the passages between these dynamics live only for 
a finite time t. On the other hand the dynamics given by (pt and g are conjugate if 
t < 00, i.e. for finite values, which can be interpreted as the stability of the dynamics. 

In Tropical geometry scaling limits result in very different objects (dynamics) which 
are far from continuous processes. Similar phenomena occur for instanton moduli 
spaces, namely compactification implies that at infinity of moduli spaces, different 
spaces appear by bubbling-off instantons. At any large scale the Chern classes are 
the same, but in the limit different Chern classes appear. 

(2) The Lotka- Volterra cell automaton (LVCA) is given by the equation: 

= + max{L, J - max{L, m^+J. 

It is shown in [10] that there is a piecewise linear map : — )■ M such that the 
corresponding dynamics $(</?) : — )■ projects to the LVCA. The corresponding 
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complex dynamics is given by: 



fL _|_ s 



For t = oo, the dynamics collapses and is described by the trivial equation 2;^^^ 



3 Intersection theory on the moduh spaces of curves 
and the matrix Airy function: 

We shall now give a mathematical application of our general considerations, namely 
to the Kontsevich-Witten model (KW) [12l[20], which originated in string theory. 

KdV flows appear both in the physical and in the biological context. One of these 
arises from scale transformations of an automaton passing through Tropical geometry 
and the other shall be Witten's generating function, i.e. the tau-function, as it was 
shown by Kontsevich. Thus from our comparisons between physical and biological 
structures, it follows naturally that the KW-model should have characteristics of an 
automaton. 

So, let us give a commutative diagram for the Kontsevich model along the lines of 
our automation program. His approach to the Witten conjecture was based on three 
key objects: 



1. One is to discretise the moduli space of Riemann surfaces via Strebel differentials 
and the thereto associated ribbon graphs. 

2. The second is to obtain the 'main identity' by use of Feynman diagram tech- 
niques. 

3. The third is to interpret Witten's generating function via random matrices by 
use of the Airy function. 



Now let or MMM stand for Miller-Morita-Mumford class and let us denote discrete 
by 'd' and ultra-discrete' by 'ud' symbohcally, such that implies a discrete r- 
function etc. . Then our program of automation of the Kontsevich model yields the 
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following commutative diagram: 

Automaton discrete dynamics continuous dynamics 







Td 




r = exp(F) 


t 




t 




t 


Airy automaton 




discrete Airy 




Airy function 


t 
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(M^ class) „d 




(M^ class) d 






t 




t 




t 










9Jl 



Let us remark that the discrete tau function can be obtained by the Miwa transform 
and that the middle objects on the vertical line have in part already been obtained 
by Chekhov ([3]). 

In [6, 7j we study another portion of the above diagram in detail, namely the one 
from the discrete MMM classes (M^ class);^ to Witten's generating function F. To do 
so, we shall combine three mathematical fields which have been developed indepen- 
dently, namely discrete integrable systems, moduh theory of Riemann surfaces, and 
the theory of mapping class groups. 

The first two are connected by passing through discrete surface theory [2|, whereas the 
last two through the group cohomology of the mapping class group. So, let us denote 
the mapping class group of genus g Riemann surfaces, which are finitely generated 
infinite groups, by F. Also, as it is known, the rational cohomology of the moduli 
space of Riemann surfaces is isomorphic to the group cohomology of F, i.e., H*{DJlg : 
Q) = H*{r : Q). Now, passing through these isomorphisms, the Poincare dual of 
the Kontsevich cycles, a* can be represented by elements of the group cohomology. 
Further, H*{r : Q) contains characteristic classes of surface bundles which are 

called Miller-Morita-Mumford classes, with known explicit representation theoretic 
constructions of the cocycles for classes. Igusa [9j obtained the leading terms of 
the Kontsevich cycles 

Now, our aim is to obtain a discrete dynamics, given by finite data, for which the 
dynamical scaling limit is Witten's generating function F. Another important point 
for us is the construction of an automatic structure {A, W, MgJ on the mapping class 
group, where A = {gi, ■ ■ ■ , gm} is a set of generators and W is an automaton with 
respect to A ( [H]). Also, the discrete Miller-Morita-Mumford class {A,W, Mg^,ip) is 
equipped with an automatic structure and a map ip : A"^^ — )• End(Z'^), such that 
these sets should admit scaling limits which permit to obtain parametrised group cocy- 
cles : F2"+^ Q and its Taylor expansion F = exp(S„e^) = Y.(k) /^fco.fei,... H^o • 
In the special case of genus g = 0, the mapping class groups are isomorphic to 
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the Braid groups divided by their centres. In relation with integrable systems, 
Zabrodin [21] found actions of the /^-matrices, which will lead in our case to the 
construction of the ip above. 

Comments: (1) On 97t^°™'^' — )■ (M^ class)^. For the construction of an automatic 
structure on the mapping class group, Mosher [15] used triangulations on punctured 
Riemann surfaces, which will be related to another arrow in our previous diagram. For 
surfaces of constant negative curvature, the angle function of asymptotic directions 
satisfies the sine-Gordon equation. On the other hand, according to the Thurston- 
Mumford approach to discretise the moduli spaces of Riemann surfaces, i.e., to give a 
cell decomposition, Kontsevich used Strebel differentials. So, it would be interesting 
to consider the relation between the angle and the differentials, from the point of view 
of Schwarzian derivatives . 

(2) The moduli spaces of Riemann surfaces 50tg,„ posses a discrete version (£OTc,,„)disc- 
Further, it is known that there is a symplectic structure u on ^g^n- Therefore, one can 
associate the symplectic Floer homology HF(iHg^„, a;), which would give an idea how 
to construct a discrete version of symplectic Floer homology directly from (^Otg,„)disc- 

(3) Yet another domain which naturally fits into our framework is quantum cosmology 
and the physics of black holes, cf. [17j 
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